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ABSTRACT 

We make use of the largest and most homogeneous sample of white dwarf/M dwarf 
(WD/dM) binaries from the Sloan Digital Sky Survey (SDSS DR7) to investigate rela- 
tions between magnetic activity, rotation, magnetic braking and age in M stars. These 
relations are studied separately for close WD/dM binaries that underwent a common 
envelope phase and thus contain tidally locked and hence rapidly rotating M dwarfs, 
and for wide WD/dM binaries that never interacted. For the wide WD/dM binary 
sample we find that the M dwarf activity fractions are significantly higher than those 
measured in single M stars of spectral type MO to M5. We attribute this effect as a con- 
sequence of M dwarfs in wide SDSS WD/dM binaries being, on average, significantly 
younger and hence more active than the field M dwarf population. The measured M 
dwarf activity fractions in wide WD/dM binaries show as well a significant increase 
from spectral types M3 to M5, where these low-mass stars become fully convective. 
This provides additional observational evidence for magnetic braking being less effi- 
cient below the fully convective boundary, in agreement with the hypothesis of fully 
convective stars having considerably longer activity lifetimes than partially convective 
stars. The M dwarfs in all our close binaries are active, independently of the spectral 
type, giving robust observational evidence for magnetic activity being enhanced due to 
fast rotation. The rotational velocities of the M dwarfs in our close binary sample are 
significantly higher than seen among field M dwarfs, however the strength of magnetic 
activity remains saturated at logLna/Lboi ~ —3.5. This unambiguously confirms the 
M dwarf saturation-type rotation-activity relation. 

Key words: Binaries: spectroscopic - stars:low-mass - stars: white dwarfs - binaries: 
close - stars: post-AGB 



1 INTRODUCTION 

The mechanisms that induce magnetic activity in field stars 
have been intensively studied and discussed over the last 
decades and this effort is now converging in a fairly consis- 
tent scenario that explains both magnetic field generation 
and angular momentum evolution in both partially and fully 
convective stars. 

Magnetic fields in partially convective (F,G,K and early 
M) stars are believed to be generated at the transition region 
between the radiative interior and the differentially rotating 
outer convective zone (tachocline) through a solar-type af2 
dynamo. In this scenario magnetic fields are generated as 

a result of a combination of differential rotation (SI effect) 

i 1 1 1 | f i 

and c yclonic convection (a effect) (Parker 1955; Lcighton 

1969 1; [Spiegel fc Zahnl Il992l ; ICharbonneaul I2OO5I ; Browningl 

2008). Lower mass M dwarfs of spectral type later than 



~M3 are fully convective and therefore do not posses a 
tacho cline, however displ ay clear signs of magnetic activity 
(e.g. iDelfosse et al.lll998h. Such stars a re suggested to ro- 
tate as rigid bodies I Barnes et al.|[2005l ). implying that the 
aQ, dynamo cannot be responsible for producing magnetic 
fields. Instead, an a 2 dynamo has been prop osed for mag- 
netic field generation in fully c onvective stars (jRaedler et al.l 
ll99d ; IChabrier fc Kiiker1l2006D . 

The efficiency of both the afi dynamo and the a 2 dy- 
namo ar e expected to be str ongly correlated with the Rossby 
number l|Noves et al.lll984l . Ro = Prat /to, with P rot the ro- 
tational period and ro the convective overturn timescale). In 
partially convective stars magnetic activity rises for decreas- 
ing Rossby numbers, implying that larger rotational veloci- 
ties lead to higher levels o f magnetic activity (for a given ro, 
lHartmann fc Noveslll987l ). Similarly, a correlation between 
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activity and rotation is also expected in fully convective stars 
in which magnetic fields may be generated by the a 2 dynamo 
l|Durnev fc Stenfldll972t IChabrier fc Kiikerl l2006h . Several 
early observational studies have explore d the connect i on be- 
tween magnetic activi ty and rotation dWilsonl 1 19661 ; iKrafj 
119671 ). More recently, iPizzolato et all (|2003t ) showed that 
the X-ray emission increases with rotation for slowly rotat- 
ing stars but saturates below a threshold of Ro — 0.1. The 
same relation seems to apply for the magnetic field strength 
ijReiners et al] |2009). Unfortunately, the phenomenon of sat- 
uration is not fully understood yet, one possible explanation 
is that sat uration occurs due to a change of dynamo con- 
figuration (|Vilhulll984l ; IWright et alj l201lh . An additional 
explanation is that saturation is reached when a certain 
maximum fraction of the available energy flux (from convec- 
tion) is converted into magnetic energy (|Christensen et al] 
2009). For ultracool dwarfs of spectral type M7-9.5 this clear 
rotati on-activity relation, however, seems to be n o longer 
valid (|Mohantv et all |2002| ; IReiners fc Basril |2010| ). More- 
over, it has to be stressed that rotation and magnetic ac- 
tivity i n single (late-type) M dwarfs might not always be 
linked l|West fc Basrill2009l ). 

Rotatio n and magnet i c act i vity decrease in time 
dSkumanicbl Il972l; iMestell Il984l; iMestel fc SpruitJ 1 19871 ; 



2008) 



as a con- 



lKawalerlll988l ; ISills et all 12000] ; IWest et al 
sequence of magnetic braking, i.e. angular momentum is 
extracted from the convective envelope and lost through 
a magnetized wind. The braking timescales stron gly de- 
pend on the spectral type (or mass) of the star (Barnes] 
120031 : n arnes fc Kim! l20ld ). which results in lower-mass 
partially convective stars spinning-down slower than their 
higher mass counterparts. Fully con vective stars have even 
longer rotational braking timescales IIReiners fc Basril 120081 ; 
iBrowning et al.ll20ld ; ISchreiber et alj|201Ch . Unfortunately 
the reason for magnetic braking being less efficient in fully 
convective stars is not completely understood. This might be 
related to a change in the magnetic field topology, which has 
been observed to switch from less ordered field structures to 
dipoles at the fully convect ive boundary going f rom partially 
to fu l ly convective stars jDonati et alj 20081; Morin et al] 
120081 ; IReiners fc Basril 120091 ; iMorin et alj l20ld) . However 
such a switch should lead to increased magnetic braking for 
fully convective stars, exactly opposite to what is indicated 
by the observed rotation rates. A new s cenario for magnetic 
brakin g has been recently developed bv lReiners fc Mohantvl 
|2012l ) who suggest that the steep transition of the braking 
rates at the fully convective boundary arises due to the sig- 
nificant drop in radius, without invoking any change in the 
dynamo theory. 

M dwarf companions to white dwarfs that form part 
of detached binaries (hereafter WD/dM binaries) offer an 
ideal test bed for studying magnetic activity, rotation and 
magnetic braking across the fully convective boundary for 
three reasons. (1) WD/dM binaries can easily be sepa- 
rated into wide binaries with orbital periods typically ex- 
ceeding 100 days, and close binaries that evolved through 
comm on envelope evolution (|Webbinklll984l ; IZorotovic et alj 
|2010|) with final orbital periods generally below one day 
( Nebot Gomez-Moran et al]|201l] ). While the M dwarfs in 
wide binaries should be unaffected by the presence of the 
white dwarf, those in close post-common envelope binaries 
(PCEBs) are tidally locked, i.e. they are rapidly rotating. 



Comparing the two samples can provide clear constraints 
on the impact of fast rotation on activity. (2) Using PCEBs 
with measured orbital periods one can furthermore directly 
relate the strength of activity to the rotational velocity of 
the M dwarf (the orbital period of the binary equals to the 
rotational period of the M dwarf which then gives the rota- 
tional velocity by adopting a spectral type-radius relation). 
(3) The age of the M dwarfs in wide WD/dM systems can 
be estimated relatively accurately using cooling tracks plus 
an initial-to-final mass relation for their white dwarf com- 
panions. This allows to test activity lifetimes that are pre- 
sumably related to magnetic braking models. 

The first studies of magnetic activity in M dwarfs that 
form part of W D/dM binaries from the Sloan Digit al Sky 
Survey (SPSS lYork et al] l200d: lAihara et alj 1201 ll ) were 
performed by ISilvestri et al] (|200d ) and more recently by 
iMorgan et all |2012l ). However, these two studies consider 
all WD/dM binaries in SDSS to b e part of close bina- 
ries, which is clearly not the case (|Schreiber et al] l20ld ; 
iRebassa-Mansergas et ahll201ll ). and their results should be 
interpreted with extreme caution. Only our dedicated ra- 
dial velocity survey of WD/dM binaries from SDSS allows 
to accurately separate wide binaries (in which the stel- 
lar components have evolved as if they were single) from 
close binaries that underwent a common envelope evolution 
and therefore contain tidally locked fast rotating M dwarfs 
|Schreiber et al.ll20ld ; lRebassa-Mansergas et alj|201ll ). This 
is fundamental to explore the relation between activity and 
rotation. The largest and most homogeneous sample of 
WD/dM binaries from SDSS inc luding radial velocity in- 
format ion has been presented in IRebassa-Mansergas et al] 
(2012aj). This sample forms the input for the present paper 
where we analyse the magnetic activity of M dwarf compo- 
nents in WD/dM binaries, thereby testing possible relations 
between age, activity and rotation. 



2 THE SDSS WD/DM BINARY CATALOGUE 

The catalogue of IRebassa-Mansergas et al] l|2012al ) consists 
of 2248 white dwarf-main sequence binaries from the data 
release 7 (DR 7) of SDSS, of which 2019 contain an M dwarf 
(WD/dM binaries). The SDSS WD/dM binary sample is de- 
scribed in more detail in this Section. We begin introducing 
our radial velocity survey that allows to separate between 
close and wide binaries. We then explain how SDSS WD/dM 
binary spectra can be used for obtaining reliable stellar pa- 
rameters of both components and finally give indications of 
how we estimate the ages of our systems. 

2.1 Close and wide WD/dM binaries 

The population of WD/dM binaries contains close binaries 
(post-common envelope binar ies or PCEBs) that underwent 
common envelope evolutio n (|Webbinkl 1 19841 ; llben fc Livid 
ll993l ; IZorotovic et al]|201ol ). and wide binaries that never in- 
teracted. The M dwarfs that form part of PCEBs are tidally 
locked and are therefore fast rotators. The properties of the 
M dwarfs that are part of wide binaries should resemble 
those of single M dwarfs as the stellar components in wide 
WD/dM binaries have evolved as if they were single stars 
jWillems fc Koldl2004| y 
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To characterize large samples of both close PCEBs and 
wide WD/dM binaries we are performing a large radial ve- 
locity follow-up program of SDSS WD/dM binaries. Our ob- 
serving strategy follows a two-step procedure. In a first step 
we obtain at least two spectra per target, separated by at 
least one night, and identify systems as PCEBs if the radial 
velocities we measure from these spectra show more than 
3cr radial velocity variation. If we do not detect significant 
radial velocity variation, we consider the system as a wide 
WD/dM binary candidate. This technique is obviously less 
sensitive towards the detection of long orbital period PCEBs 
and/or low-inclination systems and we predict that, on av- 
erage, ~16 per cent of the wide WD/dM binary candidates 
are in fact unrecognised PCEBs (|Nebot Gomez-Moran et al.l 
l201lT ). In a second step, the orbital periods of the PCEBs 
identified in this way are measured from higher cadence ra- 
dial velocity follow-up. 

So far we have identified 19 1 PCEBs and 1055 wide 
WD/ dM binary candidates (e.g. IRebassa-Mansergas et ail 
l201ll) and have measured t he orbital periods of 81 PCEBs 



e.g. Schreiber et al. 20081; [Rebassa-Mansergas et al. 2008; 



Nebot Gomez-Moran et ail l201ll : IRebassa-Mansergas et all 
2012bl ). These samples of wide and close WD/dM binaries 



offer the potential to e.g. directly compare activity fractions 
of slowly and fast rotating M dwarfs as a function of spec- 
tral type. Furthermore, the rotational velocities of the M 
dwarfs in PCEBs with known orbital periods can also be 
known (assuming a spectral type-radius relation). These M 
dwarfs are very fast rotators and should populate the sat- 
urated regime, therefore allowing us to test the saturation- 
type rotation-activity relation. 



2.2 Stellar parameters 

Reliable stellar parameters of the WD/dM binaries are 
obtained from the SDSS spectra using our decomposi- 
tion/fitting routine (|Rebassa-Mansergas et all [2007). The 
SDSS spectrum is initially fitted with a two-component 
model using a grid of observed M-dwarf templates and a grid 
of observed white dwarf templates (see Figure[T}, and the 
spectral type of the M dwarf is determined. Then, the best- 
fit M dwarf template, scaled by the appropriate flux scaling 
factor, is subtracted and the residual white dwarf s pectrum 
is fit ted with a model grid of DA white dwarfs l|Koesterl 
|2010| ) to determine the white dwarf effective temperature 
and surface gravity. From a n empirical spectral type-radhis - 
mass relation for M-dwarfs l|Rebassa-Mansergas et al.ll2007l) 
and a mass-radius relation for white dwarfs l|Bergeron et al.l 
Il995l : iFontaine et"ai1l200ll ) the masses and radii of the M 
dwarf and the white dwarf components are calculated. 

In the context of the present paper the spectral types of 
the M-dwarfs are of fundamental importance. We therefore 
compared the spectral t ypes of our t empla te fitting routine 
with those obtained by We st et al.l (120111). Fittin g several 
hundred M dwarfs from the list of West et all (|201ll ). equally 
distributed in spectral type, with our M dwarf templates and 
adopting an uncertainty of half a spectral type sub-class, we 
find tha t 98 per cent of our spectral classifications agree with 
those of IWest et al.l (|201ll ) within one spectral subclass. 
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Figure 1. SDSS spectrum of SDSSJ042200.51+073358.9, to- 
gether with the best fit white dwarf model (blue dashed line) 
and best-fit M dwarf template (red solid dotted line). 



2.3 Ages 

For a given mass, the magnetic activity of M dwarfs is 
clos ely related to their rotatio nal velocity (see for exam- 
ple ICardini fc Cassatellal [2003). If other effects, such as a 
close companion, do not affect the star, magnetic braking is 
slowing down the rotation with time. The magnetic activ- 
ity in our sample of M dwarfs that form part of wide SDSS 
WD/dM binaries should therefore primarily depend on the 
age of the system. Fortunately, for a given wide WD/dM 
binary one can estimate the age of the system, and con- 
sequently the age of the M dwarf companion, by summing 
the white dwarf cooling age to its main sequence progen- 
itor lifetimfl SPSS WD/dM binary ag es were previously 
calculated by IRebassa-Mansergas et all l|2012al ) for a sub- 
sample of SDSS WD /dM binaries. We estimate the ages for 
the wide SDSS WD/dM binary sample studied in this work 
(Section[3)l following the same strategy. 

The white dwarf co oling ages were calcula ted using the 
evolutionary tracks of of lBergeron et al.l (1 1995T PI. The initial 
main sequence progenitor mas ses were calculated using the 
initial-to-final mass r elation by Catalan et alj (|2008h and we 
used the equations of iTuffs et alj (|2004l . see their Appendix) 
to calculate their main sequence lifetime. In this exercise 
we excluded wide WD/dM binaries containing white dwarfs 
cooler than 12000 K, as the white dwarf spectral fitting over- 
estimates the mass in such systems (e.g. iKoester et al.ll2009l ; 
iTremblav et al.1 l201ll ) , which affects the resulting cooling 
age. For the same reason we only determined the age of wide 
WD/dM binaries if the error of the mass of the white dwarf is 
smaller than 0.2Mq, and given that the initial-to-final mass 
relation is only well defined for systems with white dwarfs 
masses >O.55M0, we only considered wide WD/dM binaries 
containing white dwarfs with masses above this value. 



1 Note that only reliable ages can be estimated for wide WD / dM 
binaries as PCEBs have undergone a common envelope phase and 
an initial-to-final mass relation for the white dwarf is not valid in 
these cases. 



2 We used Bergeron's updated grids available 
http://www.astro.umontreal.ca/~bergeron/CoolingModels/ 
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Figure 2. Left: the residual M dwarf spectrum of SDSSJ042200. 51+073358 (the complete spectrum can be seen in Figure[TJ after the 
subtraction of the best- fit white dwarf model (top left), the parabolic fit to the continuum around Ha (red solid line, top right) and the 
normalised spectrum (bottom). The vertical dashed lines indicate, from left to right, Ha-7A, Ha and Ha+7A where Ha=6562.76 A and 
Ha-7A < A < Ha+7A is the wavelength range in which the equivalent width of Ha is calculated. In this case, Ha emission is clearly 
detected, and we catalogue the system as active. Right: the same but for SDSSJ091333. 06+124429. 7, an inactive system in our sample. 



3 IDENTIFICATION OF ACTIVE M DWARFS 

In the context of this paper it is necessary to unambiguously 
identify both active and inactive M dwarfs among the SDSS 
WD/dM binary catalogue described in Section^ Two fea- 
tures have been predominantly used as activity indicators 
in optical spectroscopy, namely the Ca n H and K resonance 
lines and Ha emission. Ha emission is more accessible in 
WD/dM binaries because the Can H and K lines fall at 
wavelengths generally dominated by the flux of the white 
dwarf. We therefore consider a given M dwarf component 
as active if Ha emission is detected in its SDSS spectrum. 
We developed an automatic routine for this purpose that is 
described below. 

In a first step the best white dwarf model fit (see Sec- 
tionROl Figured]! is subtracted from the SDSS WD/dM bi- 
nary spectrum. In order to avoid our results being domi- 
nated by low quality data, only residual M dwarf spectra 
with a signal-to-noise ratio (S/N) of at least > 10 are usec£3, 
which results in a sample of 739 objects. For these 739 sys- 
tems the continuum flux of the residual M dwarf around Ha 
(6420A < A < 6700A) is fitted with a parabola (excluding 
the range Ha-7A < A < Ha+7A, where Ha= 6562. 76A). The 
parabolic fit is then used to normalize the spectrum. Exam- 
ples are shown in Figure[2] Once the spectrum is normalized, 
the equivalent width of Ha (EWho) is calculated within the 
range Ha-7A< A < Ha+7A ( emission in our case is indi- 
cated by a negative EW) . Active M dwarfs are then selected 
applying the following criteria: 

EWho < -0.75A (1) 
| EWhq |> 3x I eEWija | (2) 

3 Originally a S/N of at least >3 was used, however the activity 
fractions (see details in Section[4]| we obtained for late-type M 
dwarfs where very much affected by low-quality data, and we 
therefore decided to increase the S/N threshold to a least >10. 



h > 3 X N cont (3) 

where eEWHa is the EWhq error, h is the peak of the Ha 
emission above the continuum level, and N cont is the noise 
at continuum level. An M dwarf is considered inactive if, 
conversely, any of the following conditions apply: 

EWho > -0.75 A (4) 
| EWh q K 3x I cEWho I (5) 
hax N cont (6) 

The above procedure results in 373 active and 366 inac- 
tive M dwarfs among the 739 SDSS WD/dM binaries that 
form our final sample (note that for objects with multi- 
ple SDSS exposures we consider as active those systems for 
which Ha emission is detected in at least one of the SDSS 
available M dwarf residual spectra). Visual inspection of the 
active and inactive samples confirms the efficiency of our 
method, with only three per cent of miss-classifications. Note 
that these miss-classifications are due to the fact that we 
are subtracting the white dwarf contribution from the SDSS 
spectra, which may introduce uncertainties in the residual 
flux of the M dwarf if the white dwarf dominates the spectral 
energy distribution of the SDSS WDMS binary. Examples 
of active and inactive M dwarf residual spectra are shown 
in the bottom panels of Figure[2] 

From our radial velocity survey (Section l2.1[l we identify 
63 PCEBs and 599 wide WD/dM binaries among the 739 
systems. Ages are available for 118 of the 599 wide binaries 
(Sectionl2~3l>. 



4 ACTIVITY FRACTIONS 

SP SS single M dwarf activi t y frac tions have been published 
bv IWest et ail l|2004l . 120081 . l201lh for different SDSS data 
releases (DR4, DR5 and DR7 respectively), as well for in- 
dependent (non-SDSS) M dwarf samples (e.g. lHawlev et al.l 
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Figure 3. Activity fractions (black solid dots) of M dwarfs in the complete SDSS WD/dM binary sample studied in this work (left), 
in close PCEBs (middle), and in wide-separated WD/dM binaries (right). For comparison we show the single SDSS M dwarf activity 
fractions (blue empty dots). The total number of systems per spectral type bin are given at the top of each panel. Note that the total 
number of systems on the top left panel (complete sample) are not the sum of the total number of systems on the top middle and right 
panels (PCEB and wide WD/dM sub-samples respectively) at a given spectral type because we have radial velocity information only 
for a fraction of the complete sample. Note also that the error in the activity fractions increases towards later spectral types due to the 
decrease in the number of WD/dM binaries containing late type M dwarfs. 



1 19961 : [R cmcrs et al.ll2012h . All show the same general trend, 
i.e. a significant increase of the fraction of active systems 
going from partially to fully convective stars. 

The M dwarf activity fractions of our 739 SDSS 
WD/dM binaries as a function of spectral type are shown 
on the left panel of Figure[3] ( black dots) togethe r with the 
activity fractions obtained by IWest et al.l l|201ll ) for single 
M dwarfs from the SDSS. The two distributions differ signif- 
icantly at early-to-mid (MO-4) spectral types. At first glance 
this difference may not seem entirely surprising, as a signif- 
icant fraction of WD/dM binaries are in fact close PCEBs 
containing rapidly rotating M dwarfs and, as outlined in the 
introduction, fast rotation induc es magnetic activity. One 
might therefore spec ulate, as e.g. ISilvestri et alJ (|2006h and 
iMorgan et ail (|2012i ). that the increased activity fraction 
of WD/dM binaries containing early-to-mid spectral type 
M dwarfs is caused by a fraction of short orbital period 
PCEBs. However, we have shown that very few WD/dM bi- 
naries containing earl y type M dwarfs (<M3) are PCEBs 
l|Schreiber et al.M20ld ). i.e. the majority of these PCEBs 
have evolved towards shorter orbital periods due to an- 
gular momentum loss via magnetic braking and entered a 
semi-detached configurat ion becoming cataclysmic variables 
jPolitano fc WeileJl200^ 1. Therefore the fraction of PCEBs 
misclassified as wide systems in our radial velocity survey 
should be low in this spectral type range too. 

Clearly, the observed activity fraction does not only de- 
pend on the rotational velocities but also on the intrinsic 
ages of the M dwarfs. Therefore, to understand the mea- 
sured activity fractions requires to separately investigate 
close PCEBs and wide WD/dM binaries, and to take into 
account the ages of the latter population. 

4.1 PCEBs 

The S DSS PCEB orbital period distribu tion peaks at ~ 8 
hours (jNebot Gomez-Moran et al.l l201ll ). and the longest 
SPSS PCEB orbital period is slightly below 10 days 
|Rebassa-Mansergas et aTll2012bh . C onsequently M dwarfs 
in SDSS PCEBs are tidally locked (|Gladman et alJll99d ) 
and are expected to be magnetically active. 

The middle panel of Figure[3] shows that all M dwarfs 



within SDSS PCEBs are active independently of the spec- 
tral type. One might be inclined to interpret this result as 
an obvious consequence of rotation causing activity, never- 
theless it is important to keep in mind that Ha emission 
in PCEBs may arise from irradiation effects if the white 
dwarf is relatively h ot and/or the PCEB orbital period is 
relatively short (e.g. iTappert et al]|2011 a b). Therefore, we 
investigate below the effect of irradiation present among the 
63 PCEBs in our sample. For 47 of these PCEBs the or- 
bital periods are well kn own and the stellar parameters re l- 
atively well constrained dNebot Gomez-Moran et al.ll20lil ). 
With this information at hand we can easily estimate the 
influence of irradiation by calculating the ratio between the 
intrinsic flux of the M dwarf (Fi n t) and the irradiating flux 
from the white dwarf on the M dwarf's surface (Fi rr ), 



Fint _ / TdM \ 



\Twd ) 



Rwd 



(7) 



where TdM and Twd are the M dwarf and white dwarf ef- 
fective temperatures respectively, Rwu is the white dwarf 
radius, and a is the orbital separation. 

The resulting Fint /Fi rr ratios are shown in Figure|4] as 
a function of white dwarf effective temperature (top panel) 
and orbital period (bottom panel). Among the 47 PCEBs, 
eight were followed-up photometrically as part of our sur- 
vey (Section l2.1[) and inspection of their light curves reveals 
clear irradiation effects in only two objects (red solid tri- 
angles in Figure[4] cyan dots represent the remaining six 
systems). Based on the evidence drawn from these eight 
PCEBs we define Fint = 3 x Fi rr as the threshold separating 
systems dominated by irradiation from those that are dom- 
inated by the intrinsic emission of the M dwarf (horizontal 
gray dashed line in Figure^}. 93 per cent of our systems lie 
above this limit and are presumably unaffected by irradia- 
tion. Although the limit defined in that way just represents 
an empirically motivated estimate, we conclude that irradia- 
tion induced Ha plays a minor role for our sample of PCEBs 
and that all M dwarfs in PCEBs are magnetically active. 



© 2012 RAS, MNRAS 000, rjl] 



6 A. Rebassa-Mansergas et al. 




40 60 
Orbital Period [h] 



Figure 4. The M dwarf intrinsic flux / white dwarf irradiated 
flux ratio (F; nt /F; rr ) as a function of white dwarf effective tem- 
perature (top) and orbital period (bottom). The horizontal gray 
dashed line represents F; nt = 3x F; rr and is adopted to be the 
maximum limit for irradiation effects to be strong enough for H« 
to be in emission. Red solid triangles represent confirmed PCEBs 
showing signs of irradiation effects in their light curves, cyan solid 
dots PCEBs displaying no signs of irradiation in their light curves. 
The large errors obtained for several of our systems are a sim- 
ple consequence of propagating relatively larg e uncertainties on 
the e stimated effective temperatures ^Reb assa-Man sergas et aD 
120071) . 



4.2 Wide WD/dM binaries 

It is generally accepted that the stellar components in wide 
WD/dM binaries have evolved in the s ame way as if they 
were sine le stars jWillems fc Kolbll2004 ). Consequently the 
properties of the M dwarf companions in such binaries 
should be fairly similar to those of single M dwarfs, and 
naively one would expect similar activity fractions in both 
populations. This is clearly not the case (see the right panel 
of Figure[3]). 

The larger activity fractions detected in our sample 
must be caused by the fact that these M dwarfs are part 
of binary systems , but as these binarie s are fairly wide 
(Porb > 100 davs. IWillems fc Kolbl l2004l l . interactions be- 
tween the two stars can be excluded to to be the cause of the 
increased activity fractions. As mentioned in Section l2~Tl our 
radial velocity survey suffers from observational biases and 
we expect a small percentage of the wide WD/dM binary 
sample to be formed by unrecognised PCEBs. It is therefore 
important to evaluate whether this effect can explain the 
increase of the activity fractions seen in M dwarfs that are 
part of wide WD /dM binaries as compared to those of single 
M dwarfs. For this purpose we calculate in what follows the 
PCEB contamination in our wide WD/dM binary sample 
at a given spectral type. The exact number of unrecognised 
PCEBs depends on the s pectra l type of the M dwarf com- 
ponents. ISchreiber et al.l (|2010l ) demonstrated that the ob- 



served fraction of PCEBs containing early-type M dwarfs is 
~10 per cent at M0-1 spectral types, and increases to only 
~20 per cent at M2. Combining these fractions with the 
P CE B detection probability of our radial velocity survey o f 
~0.8 (|Schreiber et alj|201fj [Rebassa-Mansergas et al.ll201ll ) 
we calculate a PCEB contamination of less than five per 
cent in the MO-2 range (where the activity fractions most 
differ). Therefore this effect can certainly not explain the 
difference between the activity fractions in the single and 
wide WD/dM binary samples (right panel of Figure|3]). 

Identifying both components in a SDSS WD/dM bi- 
nary spectrum requires that neither the white dwarf nor 
the M dwarf completely outshines the other component. 
Cold white dwarfs are too faint to be detected in binaries 
with early-type M dwarfs. This implies the white dwarf to 
have an effective temperature > 1 0000 K thus excluding the 
coolest (i.e . oldest) white dwarfs (|Rebassa-Mansergas et all 
2010)Q The mass distribution of white dwarfs in wide 
WD /dM binaries is similar to the mass distribut ion of sin- 
gle white dwarfs (IRebassa-Mansergas et al.ll201ll ) implying 
that the progenitors of such white dwarfs were mainly F 
stars that spent a rather short time on the main sequence. 
In other words, the simple requirement that we need to be 
able to detect the white dwarf in the SDSS spectrum to clas- 
sify a system as WD/dM binary introduces a significant bias 
towards young systems that is not present in the samples of 
single M dwarfs. 

In order to investigate this hypothesis further we show 
in Figure[5]the activity fractions as a function of age for 118 
wide WD/dM binaries with available estimated ages (Sec- 
tion l2.3p . Two important results are extracted from inspect- 
ing this figure. (1) As expected, M dwarfs in wide WD/dM 
binaries are in general significantly y ounger than the aver- 
age age of field M dwarfs of ~ 5 Gyrs (jWest et al.ll2008l ). On 
average, younger stars rotate faster and therefore a consid- 
erably higher fraction of M dwarfs in WD/dM binaries are 
active. (2) There is a clear decrease of the activity fractions 
for older systems. Even though the errors in the activity 
fractions are fairly broad, this is in very good agreement 
with magnetic braking slowing down the rotation of the M 
dwarfs that are part of wide WD/dM binaries with time. 

It is worth mentioning that the activity fractions mea- 
sured for the systems with available ages seem to be higher 
than those obtained for the complete wide WD/dM binary 
sample (see Figure[5] and right panel of Figure[3]). This is 
simply because we can only obtain reliable ages if the white 
dwarfs are hot enough to dominate the SDSS WD/dM bi- 
nary spectra, which reflects in an effective temperature of 
> 15000 - 20000 K, depending on the spectral type of the M 
dwarf. Thus an older population containing inactive systems 



4 Note that the decrease of the number of PCEBs contain- 
ing early type M dwarfs is due to two separate reasons. First, 
the majority of these systems have already entered a semi- 
detached confi guration and thus evolved into cataclysmic vari- 
ables JSchreiber et al.ll20ld) . Second, white dwarf companions to 
early type M dwarfs need to be relatively hot to be detected in the 
SDSS spectra, therefore the S/N of the residual SDSS M dwarf 
spectra in these objects is generally lower than 10 and are ex- 
cluded from the analysis (Section[3]l . These two effects result in 
no systems with M0-1 M dwarf components being included in the 
PCEB sample studied in this work (FigurcfS] middle panel). 
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Figure 5. Fraction of active M dwarfs that form part of wide 
WD/dM binaries as a function of age. For clarity, we bin the data 
by age. The data suggests a decrease in activity fraction with 
increasing M-dwarf age, in very good agreement with rotation 
braking with time. 
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exists in our wide WD/dM binary sample for which we are 
not able to estimate the ages, making the activity fractions 
lower for the complete wide WD /dM binary sample but still 
considerably higher than those measured for field M dwarfs. 

We conclude that the age effect described above can 
explain the observ ed difference betw een the activity frac- 
tions measured by IWest et all (|201lh for single M dwarfs 
and those obtained here for the M dwarfs in wide WD/dM 
binaries. 



5 THE ROTATION- ACTIVITY RELATION 

As outlined in the introduction, rotation and magnetic activ- 
ity are strongly correlated. The strength of activity increases 
with rotation and saturates around Ro — 0.1, the exact value 
of the satur ation threshold velocit y depending on the mass 
of the star (|Pizzolato et al.l 12003). Such a saturation-type 
rotation-activity relation has been observed for the entire 
spectral type range of M dwarfs dMohantv fc B asri 2003; 
IReiners fcBasrill2010l ; IReiners et al.ll2012F) . 

In this Section we investigate the saturation-type rela- 
tion in our sample of 47 PCEBs for which we know the or- 
bital periods (Section l4.H) . For this purpose, the strength of 
magnetic activity of the 47 systems is quantified by the ratio 
between the Ha and bolometric luminosities, LHa/Lboi- This 
luminosity ratio is obtained from the meas ured EWhh val- 
ues (s ee Section[3]) following the equations of lWest fc Hawlevl 
(2008), who obtained conversion factors between the two 
considered quantities. Compared to single field stars, M 
dwarf companions in PCEBs are tidally locked and are 
therefore expected to be very fast rotators that should pop- 
ulate the saturated regime. 

Figure[6] (top panel) shows the activity strength as a 
function of rotational velocity V ro t for the 47 M dwarf com- 
panions in PCEBs with well known orbital periods (black 
solid dots) and the single M dwarfs from the catalogue of 
IReiners et alj |2012i . red triangles). Note that the latter sam- 
ple is divided into systems with accurately measured rota- 
tional velocities (solid red triangles) and systems with upper 
limits to the rotational velocities (open red triangles), and 
that these values, contrary to ours, represent projected ro- 
tational velocities (Vrot x sin i). Inspection of the top panel 
of Figure[6] clearly reveals that the M dwarfs in PCEBs pop- 



log (Rossby) 



Figure 6. Top panel: the strength of magnetic activity (evaluated 
through the Ljja/Lbol luminosity ratio) as a function of rotational 
velocity (Vrot) of 47 M dwarf s in close PCEBs. O ur data (black 
solid dots) extend the work bv lReiners et all J2012I . red triangles) 
to faster rotators and unambiguously confirm the saturation limit 
of the strength of magnetic activity at logLjia/Lboi ~ —3.5. 
Bottom panel: the strength of magnetic activity as a function 
of the Rossby number. The increased scatter of logLno/Lboi at 
log (Ro) < — 1.8 may indicate chromospheric supersaturation 



ulate the regime of the fastest rotators in which the strength 
of magnetic activity saturates at log LHa/Lboi ~ —3.5. This 
result represents a crucial and very robust confirmation of 
the proposed rotation-activity relations in M dwarfs. To fur- 
ther illustrate this point, we show in the bottom panel of 
Figure|6] the activity strength as a function of the Rossby 
nu mber Rp. The l a tter h as been calculated as described 
m Homers fc Basril i|20i0t ). All our systems fall below the 
threshold of Ro ~ 0.1, well within the saturated regime. 

A closer inspection of the bottom panel of Figure|6] re- 
veals that whilst practically all systems with log (i? );5 — 1.8 
have -3.5< logLna/Lboi ;$ —4.0, there is substantial system- 
to-sys tem scatter for sma ller Rossby numbers. Interest- 
ingly, IStauffer et aTJ l|l997l ) found a decline in X-ray sat- 
uration for stars with R < 0.01 in the coronal emission 
of young G and K dwarfs, an effect termed "supersatura- 
tion". Whilst evide nce for coronal supe rsaturation in F,G,K 
dwarfs is growing (j Wright et al.ll201ll ). coron al supersatu- 
ratio n has not been ide ntified in M dwarfs l|james et al.l 
l200d: IJeffries et all lioill) and evidence for chromospheric 
supersatur ation has not been ide n tified neither in F,G,K no r 
M dwarfs (|Marsden et al.1 120091 ; Ijackson fc Jeffries! l201Ch . 
While our data do not show a statistically significant de- 
cline of the strength of activity for very fast rotators (top 
panel of Figure[6)l , the increased scatter of the chromospheric 
activity level for systems with very small Rossby numbers 
in our PCEBs sample is intriguing in the context of super- 
saturation. One could speculate that we see first indications 
of chromospheric supersaturation in M dwarfs. However, the 
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current data are clearly not conclusive and additional mea- 
surements of chromospheric activity of M dwarfs in fast ro- 
tating PCEBs are required to test this hypothesis. 



tion of magnetic fields in low-mass stars, and add important 
information to our understanding of the relations between 
magnetic activity, rotation, age and magnetic braking in M 
stars. 



6 SUMMARY AND CONCLUSIONS 



Angular momentum evolution, rotation and age are closely 
related in low-mass stars. Magnetic field generation strongly 
depends on t he Rossby number, theref o re the rotation, of 
a giv en star (|Hartmann fc Novesl 1 19871 ; IChabrier fc Kiikerl 
2006). Thus, magnetic activity increases wit h rotation until 
its st r ength saturates for f ast rotation rates |Pizzolato et al.l 
l2003t H emers et all l2009h . Magnetic braking slows down 
the rotation of the star and conseque ntly magnetic ac- 
tivity and rotation decrease in time l|Skumanic b| 1 19721 ; 
iBrowning et al.l l201dh . The timescale of magnetic braking 
depends on the mass of the star, i.e. in fully convective 
stars magnetic braking beco mes inefficient |Reiners fc Basrl 
2008; Sch reiber et alj l201Ch and the activity lifetimes in- 
crease (|West et alJl20o3 T 

WD/dM binaries offer an excellent opportunity to in- 
dependently investigate the relations among magnetic activ- 
ity, rotation and age in M dwarfs. This is because WD/dM 
binaries can be separated into wide binaries that never in- 
teracted, and in which the properties of the M dwarf com- 
ponents should be fairly similar to those of single stars, 
as well as close binaries that underwent common envelope 
evolution (PCEBs) in whic h the M dwarf is tidally l ocked 
and rotates very rapidly dRebassa-Mansergas et al. 20 111 ; 
ISchreiber et al.ll2010l ; iNebot Gomez-Moran et al.ll201ll ). To 
test the activity-rotation-age connection with the largest 
and most homogeneous s ample of SPSS WD/dM dwarf 
binari es currently known l|Rebassa-Mansergas et all |2010|. 
l2012ah has been the main goal of this paper. Our main re- 
sults can be summarised as follows: 

• We have demonstrated that for M dwarfs in PCEBs fast 
rotation implies activity independently of the spectral type. 
This supports the connection between rotation and activity 
as a necessary ingredient for dynamo models responsible for 
magnetic field generation in both partially and fully convec- 
tive stars l|Brownindl200g ; IChabrier fc KiikerlbOOrj ) . 

• We have shown that M dwarfs in close WD/dM 
binaries populate the saturated regime (Ro < 0.1) 
and that the strength of magnetic activity saturates at 
logLHa/Lboi ~ —3.5. This represents a crucial confirma- 
tion of the saturation-type rotation-activity relation for M 
dwarfs that ar e very fast rotators, in perfect agreement with 
the results bv lReiners fc Basril ||2010T ); iReiners et al.l (|2012T l 
for single M dwarfs. 

• We have provided additional observational evidence for 
magnetic braking becoming inefficient in fully convective 
stars i.e. the activity fractions of M dwarfs in wide WD/dM 
binaries show a significant increase at the fully convective 
boundary. We have also shown that the activity fractions of 
M dwarfs in wide WD/dM binaries decrease in time, which 
proves the spin-down of these M dwarfs due to magnetic 
braking. 

We conclude the observational findings from the SDSS 
WD/dM binary sample are consistent with the current pic- 
ture of angular momentum evolution and dynamo genera- 
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